Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 2 February 2008 



(MN M?eX style file v2.2) 



Linking halo mass to galaxy luminosity 



A. Vale^* and J. P. Ostriker^ 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CBS OH A, United Kingdom 



Accepted Received ; in original form 20 February 2004 



ABSTRACT 

In this paper we present a new, essentially empirical, model for the relation between 
the mass of a dark matter halo/subhalo and the luminosity of a galaxy hosted in 
it. To estimate this, we replace the assumption of linearity between light and mass 
fluctuations with the assumption of monotonicity between galaxy light and halo or 
subhalo mass. We are enabled to proceed with this less restrictive ansatz by the 
availability of new, very high resolution dark matter simulations and more detailed 
and comprehensive global galactic luminosity functions. 

We find that the relation between halo/subhalo mass and hosted galaxy lumi- 
nosity, is fairly well fit by a double power law. That between halo mass and group 
luminosity has a shallower slope for an intermediate mass region, and is fairly well 
fit by a two branch function, with both branches double power laws. Both relations 
asymptote to L oc at low M, while at high mass the former follows L cx M'^-^^ 
and the latter L oc M°-9. 

In addition to the mass-luminosity relation, we also derive results for the occu- 
pation number, luminosity function of cluster galaxies, group luminosity function and 
multiplicity function. Then, using a prescription for the mass function of haloes in 
under/overdense regions and some further assumptions on the form of the mass den- 
sity distribution function, we further derive results for biasing between mass and light 
and mass and galaxy number, light distribution function and the void probability 
distribution. 

Our results for the most part seem to match well with observations and previous 
expectations. We feel this is a potentially powerful way of modelling the relation 
between halo mass and galaxy luminosity, since the main inputs are readily testable 
against dark matter simulation results and galaxy surveys, and the outputs are free 
from the uncertainties of physically modelling galaxy formation. 

Key v^rords: galaxies: haloes - cosmology: theory - dark matter ~ large-scale structure 
of the universe 



1 INTRODUCTION 

In recent years, N-body numerical simulations have given us 
a good understanding of dark matter structure for standard 
cosmological scenarios, while large scale observational sur- 
veys have done the same for the distribution of galaxies. In 
this way, we are now developing a good picture of how mass 
and luminosity are distributed in the universe. However, it is 
still not well known how to connect the two pictures. While 
it is well established that dark matter haloes are the hosts of 
the observed galaxies, it is still poorly understood how the 
former are related to the latter. Further, the picture is com- 
plicated by the fact that what is usually taken as a halo in 
simulations would often host multiple galaxies, especially for 
higher masses. To analyse the issue fully, it is necessary to 
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look at the halo substructure, since each subhalo can host 
a galaxy. Establishing such a link between halo mass and 
galaxy luminosity would be important because, first of all, 
it would allow us to have a direct connection between theory 
and observation, dark matter haloes and galaxies. Further, 
it could also shed some light into the theory of galaxy for- 
mation. 

There are several ways in which this problem 
can be studied. The more direct ones involve ei- 
ther numerical simulations including gas dynamics 
JWhite. Hernauist fc Springelll200ll: lYoshikawa et alJl200ll: 
'Pea rce et al.l l200lt iNaeamine et alJ l200lt iBerlind et^ 
200d 

or semi-analytical mode l s of g alaxy formation 
(Kauffmann, I^isser^^Steinmeta ^9^_jGovemat^^tab| 
19988: iKauffmann et alj I1999alllj iBenson et alJ l2000dlbf 
Sheth fc Diaferid l200ll: ISomerville et all l200ll: 
Wechsler et all I2001I: iBenson et al.l l2003al : iBerlind et^ 
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l2003r) . but, while they explicitly give the properties of 
galaxies located in a given halo, they have the added 
difficulty that many of the mechanisms involved in galaxy 
formation are poorly understood, and difficult to compute. 
Their complexity could also mask any fundamental relations 
that might be present between halo and galaxy properties. 

More indirect approaches have also been s tudied. 
The halo occupati o n distribution (HOD) model llS g^kl 
2ml lRensonll2nnn: iRullock. Wechsler SomervillJ 



g lRensonll2nnil: iRullock. Wechsler Son 
Zherie et al.ll2002l: iRerhnd fc Weinberell2002l : [ 
200,^: iMaghocchetti fc Porcianill200,'j) is based 



Rerlind et alJ 
on the prob- 
ability P{N\M) that a halo of mass M is host to N galax- 
ies. Ry specifying the P{N\M) function, along with some 
form for the distribution of dark matter and galaxies within 
each halo, it is then possible to relate different statistical 
indicators of the dark matter and galaxy distributions, such 
as correlation functions, to each other. This fully specifies 
the bias between the galaxy an d the underlying matte r dis- 
tributions. A recent paper bv iKravtsov et al.l ll2003h . has 
done a detailed study of results from simulations and re- 
lated them to the HOD model, and has concluded that the 
form of Ps{Ns\fJ,) for the subhaloes is approximately uni- 
versal, where is the subhalo mass scaled to an appropri- 
ate minimum mass. This paper also give results for the re- 
lation between galaxy a bsolute magnitude and halo circu- 
lar velocity. Other wor k ( van den Bosch. Yane & Mo '200^ 
lYang. Mo fc van den R osch 2003; Mo ct al. 2003) has taken 
this approach one step further by studying not only the 
number of galaxies associated with each halo, but also their 
luminosity, by building the conditional luminosity function, 
${L\AI)dL. This gives the number of galaxies with luminosi- 
ties in the range L ± dL/2 contained in a halo of mass Ai. 
While this work directly relates the halo mass to the galaxy 
luminosity, it does so only to the average values, lacking 
the full statistical treatment which is analysed in the HOD 
models. 

Other authors have used a slightly different method. 
Instead of trying to specify the number of galaxies in each 
halo, they treat the halo as a whole and identify it with 
a galaxy group. Then, by comparing the group luminosity 
function with the halo mass functi on, they obtain the lu- 
minosity associated with e ach halo dPeacock fc Smithll2000l : 
|m armom fc Hudsonlliooll . and also develop ways to esti- 
mate the number of galaxies hosted in a halo, thus coming 
back partly to the P(A''|A/) estimate of the HOD models. 

In the present paper, we follow a new and conceptually 
clear approach based on one simplified and testable hypothe- 
sis: there is a one to one, monotonic correspondence between 
halo/subhalo mass and resident galaxy luminosity. We might 
call this the empirical (rather than the semi-analytical) ap- 
proach because there is no attempt to physically model the 
galaxy formation process. Instead we take from observations 
the galaxy luminosity distribution and match it with the the- 
oretical halo/subhalo distribution. This has the additional 
advantage of naturally giving a lower mass threshold for 
haloes that host luminous galaxies, as the luminosity de- 
creases sharply with mass for less massive haloes. It also 
implicitly gives rise to galaxy systems, if one identifies a 
system of a massive halo and its subhaloes with the cen- 
tral galaxy and its satellites in groups and clusters. We find 
that the single assumption is very powerful and allows us 
to compute, and compare to observations many quantities 



from bias to the void distribution function to the spatial 
correlation function. 

This paper is organized as follows; in section 2 we 
present our model for the subhalo mass distribution in a 
parent halo, and build the global subhalo mass function. In 
section 3, we derive the relation between mass and luminos- 
ity, as well as some other functions such as the luminosity 
function of cluster galaxies, the group luminosity function 
and the multiplicity function. In section 4 we study how to 
apply the relation we obtain to get the light density and 
the number density of galaxies as a function of mass den- 
sity, and also obtain results for the distribution function of 
light density and the void probability function. Finally, we 
conclude in section 5. 

Throughout we have used a concordance cosmological 
model, with flm = 0.3, Qa = 0.7, h = 0.7 and ug — 0.9 
tBahcall ot aL.199a; .Soergel et alJl2003l) . 



2 SUBHALO DISTRIBUTION 

The first step in building our model is to specify the mass 
distribution of subhaloes for a given parent halo. We adopt 
the following function; 

N{m\M)dm^ A{——] exp ^ — ttt , (1) 

which gives the number of subhaloes with masses in the 
range m to m -f dm, for a parent halo of mass M. The 
normalisation A is such that the total mass in these sub- 
haloes, JJ^ m'N {m\M^dra^ is a fraction of the parent halo 
mass, x-^M (where the factor x accounts for the added mass 
of the original, unstripped, subhaloes). With this definition, 
we can write A as 



A = 



1 



/?r(2 - a) 



(2) 



This expression is motivated by recent analy- 
sis of high resolution dar k matter simulations of 



IWeUer. Ostriker fc Rodd ll2004ll . and we use for the param- 
eters the values a = 1.91, /3 = 0.39 and 7 = 0.18. Its 
re sults are simi lar to those obtained by the simulations 
of iDe Lucia et a l. ( 2004]), who find a power law fit to the 
subhalo mass function with a slope close to -2 (that is, 
in terms of equation Q a = 2). However, they also find 
that if they only include the lowest mass bins where sta- 
tistical errors are smallest, this slope is reduced to values 
around -1.9, very similar to what we use here. It further 
agrees well with the cumulative mass function derived an- 
alytically bv lOguri fc Led (|20o3) (see section below on the 
occupation numbers for further discussion), and it is also 
similar to the power law form for the subhalo number given 
in Sheth fc Jain ( 2003) , only instead of having a sharp cutoff 
at the parent halo mass, we introduce a smooth exponential 
cutoff from a mass /3M. Since it is a Schechter function, it 
is also similar to the halo mass function, and the slope a 
is close to the expected value of the slope of the halo mass 
function. 

It gives a total mass in subhaloes of 18%, close to but 
slightly higher than the values o f » 10% obtained in dif- 
feren t studies (see, for e x ample. Fib rmen. D iaferio fc Sveil 
Il998l; iGhigna et a'i]|200Ci) . iDe Lucia et all J2004l) derive a 
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lower value of around 6% at radius r2oo (where the av- 
erage overdensity is 200 times the critical), though they 
also find that this fraction can be as high as 10-15% 
in some cases. However, r2oo is smaller than the virial 
radii typically measured from the simulation results by 
IWeUer. Ostriker fc Bodel l)2004|) . which helps to explain the 
difference. 

It is worth noting that, like previo us results from sim- 
ulations and analytic a l modelling (e.g . jMgor^^^^J^^^i 
iKravtsov etal\ l2003t iDe Lucia etall 120041: lOeuri fc Led 



i Kravti 

12004^ . the shape of the subhalo mass function given by equa- 
tion Q is independent of parent halo mass. This does not 
mean, however, that one should automatically expect the 
same to be true of satellite galaxy distribution in galaxies 
like the Milky Way and clusters. In fact, as can be see from 
our results further below, the mass luminosity relation has 
very different behaviour depending on the mass of the host 
halo/subhalo; this leads to parent halos of very different 
masses (such as a cluster and a galaxy) having very differ- 
ent satellite luminosity distributions, even though the scaled 
mass functions of their subhaloes are identical. 

A very important point when analysing this expression 
is to note exactly what mass is being accounted by m. In 
fact, the distribution is valid for the present mass of the 
subhalo satellites, obtained when x = 1, which is measured 
after the tidal stripping of the outer parts of the subhalo 
in the halo potential well. However, to use this function to 
build the host distribution and to subsequently compare it 
to the galaxy luminosity function, it is necessary to use the 
original, unstripped mass of the subhalo, since only then can 
we establish a monotonic correspondence between host halo 
mass and galaxy luminosity. The factor x put into the ex- 
pression takes this mass loss into account, since we can then 
treat the mass m in equation Q as the original mass, with 
the stripped mass being then m/x with x > 1. Obviously, 
the actual factor for a given subhalo will be highly variable, 
but for reasons of simplicity and als o because it is not well 
known (see e.g. iHavashi et aljl2003l for a treatment of the 
profile of stripped subhaloes), we will be treating this as an 
average applicable to all subhaloes, and using a; = 3 in the 
present work. It should be noted that, since we are taking 
the total subhalo mass to be 18% of the parent halo mass, 
this factor cannot be more than 5, assuming that the major- 
ity of the mass in the parent halo w as built up by st r ipping 
of the subhaloes. A recent paper bv IKravtsov et al] ll2004ll 
includes a comparison between maximum and present mass 
of subhaloes in the dwarf mass range. Their results (see their 
figure 4) seem to support that the average change in mass 
does not depend on the mass of the subhalo, and an aver- 
age mass stripping factor of 3 as we adopt here seems a fair 
agreement. 

All plots and further expressions (where we will drop the 
factor x) presented refer to this original, unstripped mass; 
to invert our approximation is simply a case of dividing the 
subhalo mass by a factor of 3 to represent the actual stripped 
meiss which would be measured in a simulation. Ideally, we 
should use instead a distribution for the maximum circular 
velocity of the subhaloes. Even though this quantity also 
shows a large scatter bet ween the maximum va lue and that 
at present after stripping (IKravtsov et aljl2004ft . its relative 
change should be smaller than what would be expected for 
a given change in mass. This should be the case especially 



for subhaloes massive enough to host galaxies, which would 
have only their outer parts stripped off, and where therefore 
the inner parts which determine the peak circular velocity 
are left relatively undisturbed. This would mean that a cal- 
culation based on the cumulative number function like we 
use here would incurr a smaller error, even with a large scat- 
ter, and in later work we will transform to that variable. 

Using the expression for the subhalo number (0 to- 
gether with the halo mass function, nh{M), it is then pos- 
sible to build the global subhalo mass function, that is, the 
number density of subhaloes in a given mass range, by 

p oo 

nsh{m)= j N{m\M)nh{M)dM . (3) 

If we assume that the halo mass function has a 
Schechter form. 



n,(M)dM = c(-) exp(--) — , 

it is possible to write an analytical expression for it: 



/ m 



, _ 2C7 , „. 



Ko{2 



(4) 



(5) 



where K stands for the modified Bessel function. 

With more general mass functions, the calculation has 
to be done numerically. In the present wor k, we use the more 
accurate Sheth Tormen mass function jSheth fc TormenI 



nh{M)dM 



\ v^i IT M dM ^\ 2 ) 



(6) 



with V = \/a o(z)^(Af) ' " = O-^OT, A ^ 0.322 and q = 0.3; 
as usual, a{AI) is the variance on the mass scale M, D{z) is 
the growth factor, and Sc is the linear threshold for spherical 
collapse, which in the case of a fiat universe is Sc — 1.686. 
The result for the global subhalo mass function is shown 
in figure , where the Sheth Tormen mass function is also 
shown. 



2.1 Occupation number 

With the expression for the subhalo number, equation Q, 
it is also possible to calculate the halo occupation number, 
that is, the number of subhaloes in a halo of mass M: 



7 



(5r{2- 



rr(i 



I3M 



(7) 



Because the integral diverges, it is necessary to set a 
minimum threshold for the subhaloes. This can be set as 
the minimum mass for a halo to host a galaxy, in which 
case this occupation number corresponds to the number of 
galaxies contained in a given halo. This is one of the basic 
ingredients of HOD models, where it is taken as the average 
number of galaxies in a halo. Figure |2| shows the occupation 
number as a function of the parent halo mass, in units of 
the minimum mass considered. The most prominent feature 
is the cutoff for M/Mmin ^ 5; this is inherent to our model, 
since we consider the total mass in subhaloes to be only 18% 
of the halo mass, and we also have a cutoff in the number of 
subhaloes, given by the parameter /3. At the high mass end, 
the function is well fit by a power law, with Ns oc M'^'^^ 
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Figure 1. Mass function for haloes (solid line) and subhaloes 
(dotted line). The dashed line is the sum of the two. The parent 
halo mass function is the one proposed by Sheth Tormen (equa- 
tion|HJ, while the subhalo one is the one obtained using our model 
for the subhalo number distribution with parent halo mass. 



Figure 2. Number of subhaloes in a halo of mass M (equation 
[71 solid line). The halo mass is scaled to the minimum mass con- 
sidered for the subhaloes. The solid line shows the results for the 
model used in this paper (/3 = 0.39). The two dotted lines show 
the result of changing the subhalo cutoff (respectively for /3 = 0.2 
and /3 = 0.05). 



for (M/Mmin) ^ 1000, which is similar to what is predicted 



by most HOD models (see e.g., 


van den Bosch. Yane & Mol 


|2003| 


iMaeliocchetti & Porcianl 


20031: iKravtsov et alJl2003l). 



This behaviour is also analogous to that found in the an- 
alytically derived halo mass function of Oguri & Lcc (200^, 
who also find a slope of roughly 0.9 at high (M/Mmin), but 
a steeper value close to 1 in the lower range, similar to our 
results. The values they obtain are also similar to ours, with 
one caveat: they correspond to those we would get if we took 
the values for the stripped mass of the subhaloes, instead of 
what is plotted in figure|5| this would correspond roughly to 
diving the numbers we get by a factor of 3. Our result is also 
similar to what has been observ e d for the case of galaxies in 
a cluster, with teoch anek et al.l |2003) obtaining a relation 
for the number of galaxies with luminosity greater than L,, 
Ng{L > L*) oc MI'^ , where all quantities are normalised to 
the usual mass overdensity of 200. If we assume that their L, 
galaxy would have a mass of around 10^^ h~^MQ (see below 
for the results we obtain) , we also obtain approximately the 
same number of galaxies in a 10^^ h~^MQ parent halo. The 
overall shape of our o ccupation number is also similar to 
the results obtained bv lKravtsov et all l)2003h in their simu- 
lations, although in their case the cutoff is less pronounced, 
and in fact is only noticeable for M/Mmin ~ 1- As referred, 
however, our cutoff around M/Mmin « 5 is in fact a feature 
of our model, and something similar should be present as 
long as a cutoff in subhalo mass is introduced and there is 
an upper limit to the total mass in subhaloes. The normal- 
isation of the curves is also similar, with a value of Ns ~ 30 
for M/Mmin ~ 1000. It should however be cautioned that 
the way in which the subhalo mass is accounted may not be 
exactly the same in both cases (see discussion above). 



3 MASS LUMINOSITY RELATION 

Once we have the total mass function for haloes and their 
subhaloes, it is then possible to compare it to the galaxy 
luminosity function to obtain a mass luminosity relation. 
As noted, a few prior assumptions go into this: first, that 
the galaxies are hosted individually in the haloes or sub- 
haloes, and that each contains a single one (or, in the 
case of the parent haloes, they have one in their centre, 
plus the ones in their subhaloes); and second, that the lu- 
minosity is a monotonic function of the halo mass. The 
first assumption is supported by p revious studies of the 
group mass function. For example, iMartmez et alj ||2002|) 
obtain a group mass function from the 2dF catalogue which 
matches well with theoretical mass functions like the Sheth 
Tormen one, equation JSJ, indicating that each of their 
identified groups correspon ds to a halo. As for the sec- 
ond, Ncyrinck, Hamilton, & GnedinI i2004l) have shown that 
subhaloes identified in a set of simulations have a correla- 
tion function and power spectrum that matches the galax- 
ies in the PSCz survey, which shows it is possible to un- 
derstand the spatial distribution of galaxies by identify- 
ing galaxies brighter than a given luminosity with haloes 
larger than a certain mass. This is further justified by stud- 
ies which show that the only halo property dependent on 
the large scale environment is the mass distribution (e.g., 
iLemson fc Kauffmann|[l999^ . which coupled with the cur- 
rent understanding of galaxy formation theories should im- 
ply that we have captured most of the environmental depen- 
dence of the galaxy luminosity. 

We take the galaxy luminosity function to have the 
usual Schechter form: 
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,^(L)dL = 0,(£)%xp(-il)^. (8) 

The values of the different parameters are taken from 
the b j band 2dF galaxy luminosity function jNorberg et alJ 



1.61 X 10 



"^h^Mpc" 



Slog 10 ft 



I2002D . with $* 
— 19.66 and a — —1.21. Although this fit was determined for 
the magnitude range —16.5 > Mtj — 51og/i > —22, here we 
will be extrapolating its result to higher or lower luminosities 
as necessary. The mass luminosity relation is then calculated 
by setting the luminosity L of a galaxy hosted in a halo 
of mass M to be such that the number of galaxies with 
luminosity greater than L equals the number of haloes plus 
subhaloes with mass greater than M: 



4){L)dL ■ 



[nh{M) + n,h{M)]dM . 



(9) 



It should be noted that, because the values for the sub- 
halo mass function are generally lower than those of the halo 
one, the subhaloes make only a small contribution to this 
expression, only being important at lower mass values. This 
guarantees that a possible second generation of subhaloes 
(subhaloes of the subhaloes considered here) would not in- 
fluence this result much, since their mass function would 
have values that much lower than the original halo mass 
function, and then only at the lowest masses considered. 

It is also possible to obtain the group luminosity associ- 
ated with each halo. In order to do this, we first assume that 
a galaxy system (which would be either an isolated galaxy, 
a group or a cluster, depending on the mass of the halo con- 
taining it) can be represented by a halo and its associated 
subhaloes. The group luminosity for a halo of mass AI is 
then simply 



L{M) ■ 



L{m)N{m\M)dm. 



(10) 



There is an important point that should be considered 
when analysing this expression, and that is the possibility 
of putting a lower limit in the integral, since there may be a 
minimum mass for a halo to be host to a luminous galaxy. 
However, as can be seen in figure Q, which shows the results 
for the luminosity and group luminosity associated with a 
halo, the calculated luminosity decreases sharply with mass 
at the low end, becoming negligible for haloes with masses 
lower than approximately 10^'^h~'^MQ . Therefore, instead of 
putting in a cutoff, the value of which is not clear from what 
is known of galaxy formation, we simply use the natural drop 
off in the derived relation. 

This lower limit seems interestingly to be roughly the 
mass which is usually considered for the host ha loes of dwarf 
galaxies (for example . iThoul fc Weinberd il996l) give a lower 
limit of Vc — 30km/s to the maximum circul ar velocity of 
halo capable of hosting a dwarf galaxy; see also lStoehr et all 
(2002); Havashi et al. (2003) for their simulation results on 
the halo hosts of dwarf galaxies, where they give similar val- 
ues for the less massive of them) , but it is quite low when 
compared with the v alues used in some studies (f or example, 
in the HOD model in lBerlind fc Weinberd ll2002ll 'l. although 
this may be related to the fact that most dark matter simu- 
lations used in these studies do not reach such low masses. 
Another important aspect that can be seen in the figure is 
that the group luminosity only starts to depart from the 
halo luminosity for masses above lO^^h"^M0. This is to be 



S 10 





1 1 


1 1 ' i 


















lit 












. ^^^^ — " : 










_^// 






' '/^^^^^^^ 


j 


/ / 






/ / 
■ / / 






/ / 






- / / 






1 / 






" / / 






-/ /, 












i / 

— 1 1 . 


1 1 


1 1 . 


lO'" 10" 


10'2 10'= 


10" 10'= 




m (h-' M^) 





Figure 3. Relation between galaxy luminosity and the original 
mass of the dark matter halo which hosts it. The solid line is for 
each individual halo or subhalo, the dotted line shows the total 
group luminosity of the halo plus subhaloes system, in which case 
the X-axis m refers to the parent halo mass only. Also shown as a 
dashed line is the result when the mass is taken to be the actual 
stripped mass of the subhaloes (which is obtained by dividing 
the original mass by a factor of 3). The different points are esti- 
mates from different observational studies (see discussion in text 
for references). 



expected given the estimated mass of the host haloes of rich 
groups and clusters. This difference is the contribution of the 
subhaloes, and it is clearly seen that, although their numbers 
may not be very significant to the total mass function, they 
are very important for the luminosity of high mass haloes, 
and therefore to the luminosity function. Due to their small 
numbers when compared with subhaloes, it is to be expected 
that subsubhaloes make only a small contribution to this 
group luminosity, and that this would only be noticeable 
for the most massive parent haloes, where the subsubhaloes 
can be massive enough to be luminous. Noticeably, the group 
luminosity has a middle region of intermediate mass with a 
shallower slope, and then an upturn for high mass, which 
was also present in similar studies which used the group lu- 
minosity function (e.g.. Peacock & Smith 2000), while the 
single halo luminosity has a much shallower high end slope. 

We find that the mass luminosity relation we obtain can 
be fairly well fit by a double power law of the type: 



A- 



[m/r 



(11) 



(c+ {ni/ra'Y''y/^ ' 

For the luminosity of individual galaxies (the solid line in 
figure I^J, the parameters are A = 5.7 x 10^, m' — 10^^, 
6 = 4, c = 0.57, d = 3.72 and k = 0.23. Therefore, at the 
low mass end, we have L oc M*. In fact, if we go to very 
low luminosities, the slope will actually be slightly steeper; 
nonetheless, the region of interest essentially begins for lu- 
minosities above a few times lO^/i~'^L0 (which corresponds 
to the fainter dwarf galaxies; however, note that such low 
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luminosities are outside the range of the Schechter fit to the 
2dF luminosity hmction, so this is based on an extrapola- 
tion of the 2dF results to this luminosity range), and for 
these the slope is approximately 4. Such a slope is actually 
what you would expect from a straightforward comparison 
between the halo mass function (assuming a slope of -1.8 
at the low mass end) and the galaxy luminosity function. 
Since at the low mass end the subhaloes give an important 
contribution to the total number of hosts (see figure^, this 
is most likely a coincidence, arising from the fact that the 
total host distribution (haloes-|-subhaloes) has a slope sim- 
ilar to the halo mass function. At the high mass end, we 
essentially obtain the relation between the halo mass and 
the luminosity of the brightest cluster galaxy, which has a 
much shallower slope {L oc M^'^^). This is most likely due 
to the fact that, by construction, the mass term refers to the 
mass of the entire halo hosting the cluster, and not just to 
the mass in the region of the galaxy itself. 

The halo mass / group luminosity relation is not really 
well fit by a double power law, since it has a third region for 
middle values of mass with a shallower slope than either of 
the asymptotic values. However, we find that it is possible to 
describe it as two different double power law branches, which 
provide a fair fit to the results. Thus, for m < 10^'^/i~^Mq, 
the group luminosity is essentially the same as the luminos- 
ity of the parent halo (as these haloes do not have subhaloes 
massive enough to be luminous), so we fit it also with equa- 
tion ^2 and use the same parameters as before. For parent 
halo masses higher than 10^^/i~^Mq, we find that a simple 
double power law of the form 4.8x 10^°(m?-^-|-0.6m?''), with 
mi = m/(3.5 x lO^ft-^M©) is a good fit. Such a behaviour, 
of a curve with a relatively flat slope in the intermediate 
mass range and steeper slopes at both the low and high mass 
ends is similar to what was observed by IPeacock fc Smithl 
i200(]|) . who compared the AGS group luminosity function 
with the halo mass function, although the actual shape of 
the curve is somewhat different from what we find here. 

At the high mass end, the cluster luminosity is almost 
directly proportional to halo mass (in fact, L oc M°'^). 
This means that the resulting mass to light ratio will be al- 
most constant, rising only very slowly with halo mass, which 
matches well wit h previous re sults for the mass to light rati o 
of clusters (e.g.. iBahcall et a l. 2000; Kochanck et al.ll200.'j) . 
However, the values we obtain for the group luminosity seem 
to be smaller th an the observational resu lts. Further, the 
value derived bv lFukueita. Hogan fc Peeble s (1998) for the 
cluster blue mass to light ratio is 450 ± lQQh{M / L)q, which 
is roughly consistent with the value we obtain around 10^"^ — 
10"/i"^Mq, but smaller than what we get at IQ^^h'^MQ. 
Since the mass luminosity relation for a single halo is domi- 
nated at the high end by the parent haloes, it seems unlikely 
that this result could be much modified in the scope of our 
model, since both the halo mass function and the luminosity 
function are well known. This then means that the problem 
would lie in the subhalo distribution for these massive parent 
haloes; it would be necessary either to have more of them, 
or else for them to be more luminous. The two are in fact 
related, since an increase in the number of subhaloes causes 
an increase in the total number of hosts, and thus a decrease 
in the luminosity of a halo (especially in the low mass range, 
where subhalo number is more significant). However, our re- 
sults for the occupation number and the luminosity of lower 



mass haloes seem to be in good agreement with observations 
and prior theoretical models, which leads us to believe that 
the most likely cause of this result is the simplistic way in 
which we treated subhalo mass stripping. 

We can also compare our results with those in 
Ivan den Bosch. Yang fc Mol 1)200,?^) . who fit the mass to light 
ratio in different models they study to a double power law. 
In general, their results are in fair agreement with what we 
obtain. At the low mass end, the results are quite similar, 
with these authors obtaining a minimum of the mass to light 
ratio at a slightly lower mass. At the high mass end, they 
obtain a steeper mass to light ratio as a function of mass, 
although this is due to the flattening off we find in our re- 
sults, since in the intermediate range lO" - lO^'*/i"^M0 we 
obtain higher mass to light ratios and a steeper relation. 
This discrepancy is most likely due to the factor that these 
authors are fitting the mass to light ratio to a double power 
law, which as was discussed above does not provide a good 
fit to our results. In fact, when they adopt a different model 
with a fixed mass to light ratio at high mass, their results 
agree with ours slightly better in this intermediate region. 

Overall, our results seem a fairly good match to es- 
timates of mass taken from a range of observational re- 
sults across the entire mass range. Shown are points for the 
three most luminous of the dwarf spheroidals in the local 
group, where the luminosity was taken from the review by 
iMateo (1998i). and the mass was estimated from the results 
of lnavashi et all i2003l) (see the results in their figure 13 for 
bounds to the peak circular velocity function of unstripped 
NFW haloes estimated to be possible hosts of the local group 
dwarf spheroidals) , by assuming that the relation between 
mass and luminosity is monotonic. The dashed box repre- 
sents the relation obtained from t he weak lensing study of 
iHoekstra. Yee. fc Gladderj (|20o3|), for galaxies around L* 
(where we have taken the values applicable t o a NFW halo). 
The results for poor groups are taken from iRam ella et alJ 
I2OO2I while those for clusters are from lGirardi et aLu2000r 
where the bounds come from the two different methods the 
authors use for estimating fore/background corrections. Also 
included is a point for the Milky Way (■Cox.2000i'l . 



3.1 Luminosity function of cluster galaxies 

The major difficulty in obtaining a luminosity function for 
galaxies in clusters lies in distinguishing what haloes should 
be treated as rich clusters in our model, and which are sim- 
ply groups. We adopt the Abell definition for rich cluster, 
namely that it must have upwards of 30 objects brighter 
than ma -1-2'", where ma is the magnitude of the third bright- 
est galaxy in the cluster. Using our derived mass luminosity 
relation, it is possible to transform this magnitude thresh- 
old into a mass one. In our model, this would correspond 
to the mass of the second most massive subhalo (since we 
consider the brightest galaxy to be in the parent halo). In 
order to calculate this, we assume that its mass is given by a 
distribution which is the product of the Poisson probability 
that on average there exists a single subhalo more massive 
than it, by the probability of finding a subhalo at that mass. 
That is, its average value for a parent halo of mass Mh is 
given by 
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<MshAMh)>= / MP2(M,Mh)dM, (12) 

where P2{M, Mh) is the probabihty of the second most mas- 
sive subhalo having mass M, 

P2(M,Mh) = N(M\Mh) < N(M\Mh) > exp(- < N(M\Mh) >) , (13) " 



where N{M\Mh) is the mass distribution function of 
the subhaloes, equation Q and < N{M\Mh) >~ 

N{M'\Mh)dM' is the average number of subhaloes more 
massive than M in a parent halo of mass Mh. 

Since this mass threshold depends on the mass of the 
parent halo, in the end we obtain a condition on halo mass 
for it to be treated as a host to a rich cluster. Using our 
results, we find this to be Mn > 2.85 x lO^/i'^M©. 

If we now combine the distribution of all haloes more 
massive than this, together with their subhaloes, with the 
mass luminosity relation, we can obtain the luminosity func- 
tion of galaxies in clusters. In fact, it is possible to make 
a rough prediction for its shape at the high end by com- 
paring the luminosity of the brightest galaxy in the cluster 
(given by the solid line in figure |H1 which represents the 
luminosity of the galaxy associated with the parent halo 
itself) with the total cluster luminosity. The average lumi- 
nosity of the brightest galaxy can be estimated as Li — 

Lcf>ci{L)dL, where 4>ci{L) is the cluster galaxy luminos- 
ity function, and L' is such that on average there is only one 
galaxy more luminous than L', i .e. J^, (j>ci{L)dL = 1 (see 
lYang. Mo fc van den Boschll200 j) . With (jiciiL) normalized 
to the total cluster luminosity, Ld, it is then possible to 
extract a relation between Li and Ld, which will depend 
on the shape of the luminosity function. Taking our results 
from figureElthat Li cx A4"°-^** and Ld oc M°-^ , it is possible 
to work out that, at the high end, 4>d{L) is given by a power 
law with slope -4.21. 

We show the results we obtain in figure|l] where we have 
excluded the first brightest galaxies (in effect, we are just ac- 
counting for the galaxies in subhaloes). For c omparison, we 
include the 2dF global luminosity function of lNorberg et alJ 
i2002h . as a dotted line. We also include the luminosity func- 
tion for g alaxies in clusters in the 2dF survey, as derived by 
Ide Propris et al. ( 2003) , as the short dashed line. This is also 
given by a Schechter function, with parameters a — —1.28 
and M^j — 51ogio/i ~ —20.07. Since the normalisation is not 
given, we adjusted this to have values in the same range as 
the ones we obtain. Also shown, as the long dashed line, is 
a double power law with the same low end slope, normalisa- 
tion, and L, as the ones used to plot the previous curve, but 
where the exponential cutoff has been replaced with a power 
law with the same slope as calculated above, namely —4.21. 
Our results show a good agreement to those obtained from 
2dF, which is even better if a double power law is considered 
instead of a Schechter function. 

3.2 Group luminosity function 

Another way to check our results is to build the group lu- 
minosity function. This is the equivalent of the galaxy lu- 
minosity function, only applied to groups, and it is usually 
obtained from galaxy catalogues by building groups of grav- 
itationally bound galaxies. In our case, we start with the 
same assumptions for groups described above, and use the 
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Figure 4. Luminosity function of galaxies in rich clusters. Shown 
as the solid line is our model, where we have excluded the first 
brightest galaxies, and where we take an effective threshold of 
Mh > 2.85 X lO"/i-lM0 for the mass of a halo associated with 
a rich cluster. The dotted line shows the 2dF global luminosity 
function, the short dashed one the luminosity function of clus- 
ter galaxies also from 2dF, and the long dashed line is a double 
power with the same parameters as the previous function, but 
with 4>d{L)dL oc L-^-^i at high L. 



halo mass function and the relation between halo mass and 
group luminosity given by equation llOII so that 

MLa)dL, = nu{M{L,))^^dL, . (14) 

We show our r esult in figure El along with the AGS 
for the CfA survev llMoore. Frenk fc Whitdll993l) and the 
VSLF jMarinoni. Hudson fc Giuri cin '20021 group luminos- 
ity functions. While at higher luminosities our model would 
seem to underpredict the abundance of groups, it reproduces 
the slope of the AGS luminosity function quite well. This is 
undoubtedly related to the fact that our calculated group 
luminosity is somewhat lower than the observed values (see 
figure|21l; an increase in group luminosity would cause a shift 
to the right of our curve, providing a better match for the 
observational results. It should also be noted that the inclu- 
sion of subsubhaloes would probably tend to slightly increase 
the values at the high end, relatively to the low end, since 
most of these subsubhaloes would be dark for low mass par- 
ent haloes, but be more massive and therefore contribute a 
small luminosity to the group total in the case of massive 
parent haloes. It is also curious to note that, since the fits 
shown in the figure are a double power law in the case of the 
AGS and a Schechter function in the case of the VSLF, our 
result would probably not be well fit by a Schechter type 
function. 
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Figure 5. Group luminosity function. The solid line repre- 
sents our result using equation I14i : the two other lines are 
observational group luminosity fun ctions extracted from galaxy 
catalogues, respectively the AGS (Moore, Frcnk & White 1993i) 
(dotted) and the VSLF iMarinoni. Hudson fc Giuricia .20021) 
(dashed). 



3.3 Multiplicity function 

It is also possible to derive the multiplicity function,tlie num- 
ber of groups/clusters as a function of their richness, by a 
process similar to the one used for the group luminosity 
function. Only in this case, we use the occupation number 
(plus one to account for the central galaxy, which is hosted 
by the parent halo itself) instead of the group luminosity. 
Then, using an expression similar to 1141 with the subhalo 
number A'^s put in place of the group luminosity Lg, we ob- 
tain the result shown in figure [5] The main point to take 
into account in this calculation is, as was referred in our 
above discussion of the occupation number, the need to in- 
troduce a minimum mass for the subhaloes. In the case of 
the figure shown, this was taken to be the mass equivalent 
to Mb ~ —19.4, but in general, and to compare to results 
from observational studies, this should be set to equal the 
minimum luminosity considered for objects in the observa- 
tions. There is also a sharp upturn at A'^ = 1 due to all the 
haloes that are massive enough to be considered to host a 
galaxy, but not enough so that they can have subhaloes with 
galaxies in them. 

Figure|S|shows our results in comparison with some ob- 
se rvational data. Th e poin ts were taken from the analysis 
of IPeacock fc SmithI (|200ff) . while the two lines were con- 
structed from the two group luminosity functions shown in 
figure 13 by using the luminosity functions of the two sur- 
veys to relate total group luminosity to the number of galax- 
ies above the magnitude threshold. The galaxy numbers for 
both the points and the results derived from the group lu- 
minosity functions were then multiplied by a factor of 0.66 
to take into account the difference in radius between friends 




Figure 6. Multiplicity function derived from our model (solid 
line). This equates the number of galaxies in a group with the 
number of subhaloes in a halo, plus one to account for a central 
galaxy hosted by the halo itself. The minimum mass taken for 
the subhaloes was in this case the equivalent to a magnitude of 
Mb = —1 9.4. The different point s are taken from the results 
derived by iPeacock fc SmittJ i200(t from the CfA survey, while 
the two additional lines are derived from the group luminosity 
functions shown in figure |F] where the magnitude limit for all of 
these is also Mb = —19.4. 

of fr iends estimates and th e usual definition of virial radius 
(see iKochanek et al.ll2003h . 

Even though our results lie in the range between those 
estimated from the VSLF and the AGS group luminosity 
functions, they seem to be a bit lower than the observa- 
tional points. Since the group abundance is directly related 
to the halo mass function, which is well known, this discrep- 
ancy must be caused by the occupation number. A slightly 
higher occupation number would shift our curve further to 
the right in the figure, bringing it into better agreement 
with the observational data. Incidently, a higher occupation 
number would also bring our results for the group luminos- 
ity and the related group lumininosity function into better 
agreement with observational values, so we believe that this 
is where the problem lies. Our calculated occupation num- 
bers seem to be in fair agreement with expectations; how- 
ever, they depend sensitively on our prescription for subhalo 
stripping, where we adopted a simplistic approach. A more 
detailed and correct model for this effect would most likely 
produce better results. 



4 BIAS AND PROBABILITY FUNCTIONS 

4.1 Mass, light and number densities 

To derive a relation between the mass and light densities, 
we first need to obtain the mass function of haloes for re- 
gions of different density. In order to do this, we follow the 
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method outlined in iGottlober et all (l2003ll . In essence, we 
treat the evolution of the matter distribution in a region 
as if it derived from a universe with the same cosmological 
parameters as that region, namely the same average mass 
density as is found in the local region. We start by labelling 
each region by its average mass density, p, to which corre- 
sponds a given value of the parameter ^Im ~ p/Pc, with pc 
the critical density. Overdense regions will have ilm > ^m, 
where the barred quantities refer to the background uni- 
verse, while voids will have the opposite. 

The growth of perturbations will be affected by the dif- 
ferent matter content, and this is reflected by changing the 
normalisation of the power spectrum according to the differ- 
ential of the growth factor, D{z), relative to the background; 



as — (Tg- 



(15) 



D{0) D{z,) ' 

where Zi is some initial redshift for which the fluctuations 
were equal in the background and in this region; in the 
present case, we use Zi = 1000. We also use the usual 
normalisation for the growth factor, -D(O) = 1. To obtain 
the mass function for a region with a particular average 
density, np{m)dm, we then apply this new normalisation 
to the Sheth Tormen mass function given in equation ©, 
changing the p,„ term as appropriate, and also the value for 
^e, which has a small depende nce on Qm (see for example 
iNavarro. Frenk. fc Whit j|l997l) . 

An important point is that this prescription does not 
utilise an explicit smoothing radius for the region consid- 
ered. Instead, these are labelled by their average density. 
However, when regions of limited size are considered, there 
exists a maximum mass for objects in them, given by the 
total mass they contain. As the above prescription does 
not take this into account, we introduce an additional term 
to compensate for this effect. Therefore, after we obtain 
the mass function, we put in a further cutoff of the form 
exp[—r]{in/ms)'^], where m is the mass of the halo and 
mg = 4iT{l + Sm)pR^ /3 the total mass in the region of radius 
R and average density {1 + S,n)p. The parameter ry allows for 
some tuning of the actual cutoff, and in the present case we 
use rj — 1. Once this cutoff is introduced, it becomes neces- 
sary to renormalise the mass function, so that it still gives 
the appropriate density when the mass of all haloes is cal- 
culated. Two examples of the resultant mass functions are 
shown in figure |7| with curves for regions with the average 
density, 1/10 of the average and 10 times the average, for 
two different radii. 

Once the mass function for an over or under dense re- 
gion has been determined, it is then simple to find the equiv- 
alent light density, pL , by using the group luminosity corre- 
spondent to each halo, given by equation IIUII : 



PL 



Lg{M)nh{M)dM . 



(16) 



Our result is shown in figure |H| The results for R — 
8/!-~^Mpc and R — 4/i~^Mpc are quite similar, with the 
curve for R = l/i~^Mpc being different. This is due to the 
large suppression of high mass haloes, even at high over- 
densities, which can be seen in figure |7| which alters the 
proportion of luminous to non-luminous haloes in favour of 
the latter. T he curve we obt ain for R = 8/i~^Mpc is similar 
to the one in lMo et al.N2003h . where the authors have used a 
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Figure 7. Mass functions for regions of different size and average 
density, calculated according to the prescription described in the 
text. The upper panel shows the results for R = lh~^Mpc, the 
lower for 8/i~^Mpc. In both, the solid curves are the background 
Sheth Tormen mass function, the dotted ones are for densities 
of(top to bottom): lOp, Ip, O.lp. 



mass function for different densities derived from simulation 
results, to which they then apply the conditional luminosity 
function. The agreement is quite good with our result, even 
though we used a theoretical model for the mass function 
instead of taking it from simulations. 

The most striking feature in these results is the sharp 
decline in the light density in underdense regions. This would 
mean that the majority of void regions would be quite dark. 
Another important feature is that for the most part 5l < Sm, 
and the slope of the curve at high density is lower than 1. 
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Figure 8. Light density, 5^,, as a function of mass density, (5m- 
Both axis are scaled to the average background values, so regions 
with 1 + (5m > 1 are overdense and regions with 1 + 5m < 1 
are underdense. The different lines are for different smoothing 
lengths: R = Ih-^Mpc (solid), R = ih-'^Mpc (dotted), R = 
8/i~^Mpc (short dashed). The long dashed line marks (5^ = (5m- 
T he dashed dotted line is the fit to the same relation mentioned 
in lMo et all |20o3). 



The overall results (e.g., darkness of the voi ds) are similar 
to tho se from hydrodynamical simulations bv lOstriker et all 
i2003h . But the simulation results show a sharper cutoff at 
low mass density and hght being more overdense than mass 
for overdense regions and a slope to the relation greater than 
on e. This dif f erence is the more remarkable since the results 
of iMo et al.l ll2003t) are similar to ours while using a mass 
function derived from simulations, which otherwise might be 
considered the most likely origi n of the discrepancy. How- 
ever, previous results shown in iBahcall et al.l ll200Cf) . also 
taken from simulations, seem to be more in line with what 
we obtain in the present work. These authors in fact give an 
explanation for the apparent increasing antibias at higher 
mass densities: at low redshift these regions usually consist 
of rich clusters and superclusters, whose stellar population 
tends to be old. Therefore, young blue stars tend to be rare 
and consequently the total luminosity in the blue band is 
lower than what could be expect from their high mass, giv- 
ing rise to the slight antibias. It is therefore likely that this 
difference is coming from the way in which hght is being 
counted in the two methods. In fact, one problem with the 
approach presented here is that it breaks down for haloes 
whose size is comparable to the size of the region being con- 
sidered, since it is then possible for a sphere to encompass 
an outer region of the halo where a subhalo and a galaxy lie, 
and which is theref ore luminous but is co nsidered dark in our 
prescription. Since lOstriker et al.l ll2003h study the light dis- 
tribution directly, they would account for the light present 
in such a situation. This effect would be the more notice- 
able for high density regions and small smoothing lengths. 
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Figure 9. Same as figure (SJcurves have same label), but for 
Sg, the overdensity in the number of galaxies. The minimum 
threshold was set at halo mass corresponding to a luminosity 
L = Lf The average number of galaxies brighter than this 
limit in the background, taken from the luminosity function, is 
3-22 X lO^^/i^Mpc-^, 

and in fact the results presented in lOstriker et al.l ll2003f) 
are for R — l/i~^Mpc, while the major discrepancy with 
our results occurs for high density regions. The differences 
in underdense regions can most likely be explained by lim- 
ited resolution effects. 

We also derive results for the bias between galaxy num- 
bers and the dark matter mass density. The procedure is sim- 
ilar to the adopted for the luminosity, replacing the Lg(M) 
term in (I16II with the number of galaxies in a halo of mass M, 
i.e. (1 -I- Ns{M, nimin)), corresponding to the central galaxy 
hosted in the halo and the ones present in the subhaloes. The 
integral should also be taken with a lower limit at irimin- The 
minimum mass rrimin corresponds to the intended cutoff in 
mass (or luminosity) for the galaxies to be counted. Our re- 
sult is shown in figure |5| The general behaviour is similar 
to what was observed for the luminosity. According to these 
results, the galaxy number distribution would be a fairly 
unbiased tracer of the underlying dark matter mass in mod- 
erately overdense regions, as would b e expected from previ- 
ous studies (e.g.. iBenson et al.ll200lfl . but the galaxy num- 
bers severely underestimate the mass density in underdense 
regions and consequently assuming a direct proportionality 
can le ad to a significant underestimate of f2m iOstriker et alJ 
l2003ft . 

4.2 Distribution functions 

Using the relations obtained in the previous section, it is 
also possible to derive a distribution function for the light 
density (i.e., the probability of a region of a given size hav- 
ing a certain average light density, the light PDF) . To make 
this calculation, besides the relation between the mass and 
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light densities it is also necessary to have the mass distri- 
bution function. In general, this can be obtained from dark 
matter simulations. For the present work, we will use the 
approximation of treating this as a lognormal, with an ap- 
propriate dispersion dependent on the radius being consid- 
ered and on the cosmological model. It has been known for 
some time that such a function is a fair approximation to 
the real mass distribution function, at least in re gions of 



moderate over /under densi ty (|Cole^^^one^^^^ : see also 
iKavo. Taruva fc Sutol200ll and lostriker et ali200a for some 
recent analysis on this subject). 

The form of the lognormal is given by: 



1 



V2 



TTOj^ y 



-exp 



2aj2 



(17) 



where y = 1 + 5m and the parameter of the distribution ui 
is related to the smoothing length and the variance a of the 
mass density field by 



2 



: ln(l + ai) 



(18) 



Here, the variance a refers to the full, nonlinear spectrum. 
To relate this to the linear variance, and hence to the cos- 
mologi cal parameters we use, we follow the prescription out- 
lined in lPeacock fc Doddd il99d) to convert the linear power 
spectrum of mass density fluctuations to the nonlinear one. 

To obtain the light distribution function, we then follow 
a similar procedure to the one used for the group luminosity 
function: 

9U)dj = f{p)^dj, (19) 
dj 

where we use the notation j = 1 + Sl and g(j)dj is the 
light distribution function. Such a transformation assumes 
a monotonic dependence of the light density with the mass 
density, which is true of our results (see flgure|HJ. Another 
important point to note is that this function is not nor- 
malised to one, since some of the underdense areas (in terms 
of ma ss) are non luminous (see discussion in lOstriker et all 
|20o3). Our result is shown in figure IIUI The general be- 
haviour is what would be expected: the distribution func- 
tions are similar in the overdense regions, but in general 
differ substantially in the underdense regions, with the light 
distribution function in general not well fit by a lognormal. 



4.3 Void probability function 

The void probability function is defined as the probability of 
having no galaxies in a sphere of radius R. Although voids 
are a particularly striking feature in the galaxy distribution, 
their study is impaired by their large size and low number of 
galaxies present in the observational case, and by resolution 
difficulties when the study is conducted through numerical 
simulations including galaxy formation. Due to these factors, 
in the past studies of the voids in the galax y distribution 
have been relatively few (see, for e xample, |Vo£Blev et_al 
Il994| for observat i onal re suhs, and iMathis fc White! 12002 
and iBenson et al] l2003bl for comparisons with theoretical 
results; there are also two recent papers with observa- 
tional results for the 2dF survey, iHovle fc Voeelevl l2003t 
ICroton et al.ll20o3 ^. Despite this, the void probability is a 
powerful tool in the analysis of a particular model of galaxy 
formation, since it probes a highly non-linear regime and is 




Figure 10. Distribution functions for mass (solid) and liglit (dot- 
ted) density. The three curves correspond to the three smoothing 
lengths used previously: 1, 4 and 8 h^Mpc~^ , with peak height 
increasing with radius. The x-axis represents either mass or light 
overdensity, depending on the curve. The distribution functions 
are shown scaled appropriately for logaritlimic binning. 



not derivable from the low order correlation functions (in 
fact, it depends on all of the N-point correlation functions). 

In the case of our model, we can determine the void 
probability function by combining the dark matter distribu- 
tion function with the expected galaxy number in regions of 
a given size. Thus the number density of galaxies in a region 
of radius R and average density p, given a minimum mass 
cutoff of Mmin will be: 



wh(p) = 



{l + Ns{M))np,R{M)dM , 



(20) 



where np^fi{M) is the mass function in a region of size R and 
average density p (see discussion in section 4.1, also figure 
0, and Ns{M) is the number of subhaloes in a halo of mass 
M. The number is then obtained by multiplying this density 
by the region volume, V = 4iyR'^ /3. We treat this as the av- 
erage number an d assume a Pois s onian distribution around 
this average (see iKravtsov et alJ i2003r) . where the authors 
have shown that the full HOD is consistent with a Poisson 
distribution for large host masses), so that the probability 
of having galaxies will in fact be exp{—VnR{p)) . Finally, 
the void probability function is obtained by integrating this 
number times the probability of a region having this given 
density. 



Po{R)= / eM-VnR(p))fR{p)dp, 



(21) 



where fnip) is the probability of a region of radius R hav- 
ing density p, which we take to be given by the lognormal 
distribution, equation (1171 . 

Our results are shown in figure 1111 where we also 
show observational results for the 2dF survey taken from 
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Figure 11. Void probability function. The solid line represents 
our concordance model. The po ints are from th e combined CfA-1 
and Cf A-2 surveys (taken from lBenson et an t2003b. original re- 
sults in lVogelev et al .11994) . The observations are for a magnitude 
limit of Mg — 51og/i ^ —19.5, and the corresponding mass was 
taken as the minimum mass limit in the calculations (equation 

iHovle fc VogelevI ll2003l) . for both the NGP and the SGP. 
At low radius, we obtain a good agreement, but our model 
would seem to overestimate the probability at large radii. 
Our results are also in better agreement w ith the values ob- 
served for the SGP than for the NGP, since iHovle fc VogelevI 
i2003h have found the former to be somewhat emptier than 
the latter. There is an additional effect that may go some 
way towards explaining this discrepancy, and that is cor- 
rections for the peculiar velocity distortions. These would 
lead, especially for the larger radii, to the inclusion in voids 
of galaxies that are not actually there, due to the smearing 
effect caused by peculiar velocities to the positions of galax- 
ies in redshift space. This effect is not symmetric, since the 
voids are areas of low mass density and therefore the galaxies 
in them will have low velocity dispersions. The final result 
of this would be to overpopulate the voids and consequently 
to underestimate the measured void probability function, 
which would lead to a better agreement with our results. 



5 CONCLUSION 

In this paper we present a new model for relating halo 
mass to hosted galaxy luminosity, based on the dark mat- 
ter substructure. We feel this is a potentially powerful way 
to approach this problem, since it is based on two main 
inputs, the halo/subhalo distribution and the galaxy lumi- 
nosity function, which can be well tested and adjusted to 
results from simulations and surveys, respectively. Addition- 
ally, the model requires only one further assumption: that 
there is a one to one and monotonic relation between imbed- 
ded halo/subhalo mass and galaxy luminosity. This is more 



explicit but in general agreement with the general assump- 
tions made in past studies and general assumptions made 
in similar work of the same subject. It is a much less re- 
strictive assumption than the still sometimes utilized ansatz 
of a linear "bias" between galaxy numbers and dark matter 
density. 

We have shown how, starting with a prescription we 
describe for the subhalo mass distribution in a parent halo, it 
is possible to obtain a relation for the luminosity of a hosted 
galaxy, as well as the group luminosity when the system of a 
halo and its subhaloes is identified with groups or clusters of 
galaxies. The subsequent results appear to match well with 
general assumptions of the behaviour of such a relation, as 
well as results for the mass and luminosity at scales of dwarf 
galaxies, L* galaxies and massive clusters. 

From this model, it is then possible to derive many 
quantities that can be directly compared with further sim- 
ulation or observational results. We have shown four exam- 
ples of this, namely the occupation number, the luminosity 
function of cluster galaxies, the group luminosity function 
and the multiplicity function. Using further assumptions on 
how to calculate the mass function for regions of different 
average densities, and on the shape of the dark matter dis- 
tribution function, we have also obtained a relation between 
mass and light and number densities for different smoothing 
lengths close to what is expected from previous bias studies. 
We also obtain the distribution function of light density, and 
the void probability function. The latter is a powerful addi- 
tional test, since it probes a highly non-linear regime, and 
our results seem to match well with previous observations. 

The major difficulty with the model as it is presented 
here has to do with the identification of the mass of the sub- 
haloes. To be able to apply the monotonic correspondence 
to the galaxy luminosity, the original mass in the subhaloes 
has to be taken into account. However, the mass distribution 
we use is measured for the stripped mass of the subhaloes 
in their parent halo. In order to build the model, we took 
the approximation of taking an average for the stripping 
factor. This is, however, a not wholly satisfactory approach, 
since in general the stripping history will be highly variable 
from subhalo to subhalo, and also in different parent haloes. 
Ideally, we would like to use the maximum circular veloc- 
ity instead of the mass to identify the subhaloes, since this 
quantity should less sensitive to stripping for the relatively 
massive subhaloes which can host galaxies. Unfortunately, 
we do not have at present a good distribution for subhalo 
abundance as a function of the maximum circular veloc- 
ity, although we wish to study this further in the future. 
Nonetheless, although somewhat crude, the approximation 
we took is reasonable, and since this work deals with statis- 
tical averages for most of the quantities, it does fit somewhat 
well into the whole structure. More importantly, the results 
we obtain seem for the most part to match well with obser- 
vational data. In future work which takes into account the 
statistical variation, however, this factor will certainly be of 
great importance. 
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